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ABSTRACT 

The Relationship between the Neumann system and the Jacobi system in ar- 
bitrary dimensions is elucidated from the point of view of constrained Hamilto- 
nian systems. Dirac brackets for canonical variables of both systems are derived 
from the constrained Hamiltonians. The geodesic equations corresponding to the 
Rosochatius system are studied as an application of our method. As a consequence 
a new class of nonlinear integrable equations is derived along with their conserved 
quantities. 



1. Introduction 

To deal with Hamiltonian integrable systems, we usually rely on some basic 
procedures such as determination of action-angle variables, separation of variables, 
and integration of equations of motion. In addition, it is often very instructive to 
study relations which exist among some superficially different mechanical systems. 

For instance, it is known that there is a close relationship between the Neumann 
system [1] and the Jacobi system [2]. The Neumann system describes the motion 
of a point particle on the N — 1 sphere, S*^" 1 , under the influence of a quadratic 
potential in iV-dimensional space. Although the Neumann system has been known 
for more than 100 years as an integrable nonlinear mechanical system, it is still 
under active study [3-7]. On the other hand, the Jacobi system describes the 
geodesic flow on an ellipsoid. At first sight these two mechanical systems are 
different. However, Knorrer [7] found the mapping from the Jacobi system onto 
the Neumann system, showing that the two systems are essentially equivalent. 

A more complicated problem was recently studied by Adler and van Moerbeke 
[8]. They found that there are rational maps transforming the Kowalewski top and 
the Henon-Heiles system into Manakov's geodesic flow on SO (A). 

In the present work we deal with some classical integrable systems such as 
the Jacobi, Neumann, and Rosochatius systems laying special emphasis on the 
relationship which exists among them. In sections 2 and 3, we deal with the Jacobi 
and the Neumann systems in the framework of the constrained Hamilton formalism 
with a set of second-class constraints. We calculate the Dirac brackets for canonical 
variables in the Jacobi and the Neumann systems. We then transform variables 
so as to reduce second-class constraints to first-class ones. As a consequence, both 
systems acquire gauge freedoms. 

The relationship between the Neumann system and the Jacobi system is clari- 
fied by making use of the Gauss map introduced by Knorrer [7]. As an application 
of our method, we then study in section 4 the geodesic equations corresponding 
to the Rosochatius system [9,10], which is also known as an integrable nonlinear 
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system. Hamilton's equations for the iV-dimensional Rosochatius system can be 
obtained from the 2iV-dimensional Neumann system. By the same token, we found 
equations corresponding to the geodesic equations for the Rosochatius system from 
those of the 2iV-dimensional Jacobi system. We show that the equations obtained 
are integrable by finding integrals of motion explicitly. The last section is devoted 
to summary and concluding remarks. 

2. The Neumann system 

In this section we deal with the famous Neumann system [1] in arbitrary dimen- 
sions in the framework of the constrained Hamiltonian formalism with second-class 
constraints. Recently, various aspects of the Neumann system are investigated by 
many authors [3-7], because the Neumann system is interesting as one of the sim- 
plest integrable non-linear systems and it is closely connected to other non-linear 
mechanical systems, such as those of the KdV, the Kowalewski top etc.. Our aim 
here is to clarify the Poisson structure of the Neumann system. 

Let (x\, ...,xn, v\, vjq) be canonical coordinates, so that the Poisson brackets 
are given by {xi, Vj}p = 5ij, {xi, Xj}p = {vi, Vj}p = 0, and the symplectic form 
iv is written as iv = J^fcLi ^ x k A d v k- The Hamiltonian for the Neumann system is 



written as 



(i) 



where £ and r] are Lagrange multipliers for the second-class constraints 




k=i 



(2) 



N 




k=l 
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and the classical Hamiltonian Jijq for the Neumann system is given by 

1 N 

k=l 

with 0) and < a\ < a2 < • • • < ajv being real constant parameters. 
The Hamilton equations derived from (1) are 

ik-Vk- r]x k = 0, 

Vk + a k x k + &k + VVk = 0, 

N 

k=l 

N 

^x k v k = 0. 
k=i 

It follows from (4) and (6) that 

N N 

^2 Xk±k = S xkVk + ^ r ' 2 = °' 

k=l k=\ 

so that 

7] = 0. 

Also we have from (4), (5) and (7) 

N N 

^2(x k v k + v k x k ) = ^(vl - a k x\) - £ r 2 = 0. 
k=l k=l 

It follows from (4), (5) and (9) that 

x k = v k = -a k x k - £x k . 



Combining (10) with (11) one obtains 



N 



x k = -a k x k 2 ZJ 



•2 _ 2\ 
X j CLj Xj) 



(12) 



The Dirac brackets [12,13] evaluated from (1) are given by 



{xi, xj} D = 0, (13) 

vj}d = kj - ^y-, (14) 

X 

r t {%i v j ~~ x jVi) , ■. 

= — , (15) 

X 



where x 2 = J2k=i 



N 2 



Let us now introduce new variables m, (i — 1, ...N), defined by 



Vi = yi-^(x-y), (16) 

x z 



where x ■ y = XlfcLi ^feZ/ifc- By this change of variables, one obtains the identity 



TV 

X • V = 

k=l 



^x k v k = 0. (17) 



Therefore, the term proportional to rj in the Hamiltonian (1) vanishes identically. 
We see that (7) also turns out to be an identity. It should be noted here that the 
change of variables gives rise to a qualitative change of the class of constraints. 
Namely, although Q\ and Q2 were originally second-class constraints, after the 
change of variables f^2 vanishes identically and fii turns out to be a first-class 
constraint with the nonvanishing Lagrange-multiplier £. As a consequence, the 
system acquires gauge freedom and £ becomes the gauge parameter. 
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In order to see the gauge freedom explicitly, let us introduce the Lagrangian C 
for the Neumann system by the Legendre transform 

TV 

c = J2*kv k -n. (is) 
fe=i 

Setting Q2 = and substituting (16) into (18), we obtain 

TV TV TV 

k=l fc>/ fe=l 

where A = C — ((# • y)/x 2 ), £ — £ , = x^y/ — x/y^, and total derivative terms 
with respect to time have been deleted from (19). 

The Euler equations obtained from (19) are 

(x ' y^j 

Xk = Vk — x k , (20) 

TV 

^2 XkXk = 0, (21) 
fc=i 

_ /7, -I- - 

It follows from (20) and (22) that 



y k = -\x k - a k x k + ^(x ■ y)y k . (22) 



TV 

x k = -a k x k - ^^2(xf - aixf), (23) 
1=1 



where we have used the identity 

TV 



^4 = !/ 2_(£_|)l (24) 

fc=l 



Eq.(21) implies that 

TV 

x 2 = ^^x k = const.. (25) 
fc=i 

On the other hand, changing variables from v k to y k by (16), we have from 
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(13), (14) and (15) that 



{x h Xj} P = 0, (26) 
{xi, Vj}p = Sij, (27) 
{Vi, Vj} P = J X ^~ X ^\ (28) 

It should be noted that the Dirac brackets (13), (14), and (15) turn out to be the 
simple Poisson brackets here. 

Since is invariant under the gauge transformation, 

(xk, Vk) => ( x k, Vk + Ax fc ), (29) 
where A is any function of x and y, the Hamiltonian 

1 N 1 N 

Hn = ^2 J kl + 2 Yl akX k ( 30 ) 
k>l k=\ 

is also invariant under the gauge transformation (29), and so does the equation 
of motion (23). Although the parameter A in (19) plays a role of a Lagrange 
multiplier in dealing with the constraint Y^k=\ x k%k — imposed on the system, 
it can be regarded, on the other hand, as a gauge parameter. One notices that the 
change of variables from to (k = 1, N), gives rise to a transmutation of 
constraints from the second-class to the first-class ones, so that a gauge freedom 
appears through the gauge parameter A in the system. 

The Neumann system is Liouville integrable, because there exist N — 1 inde- 
pendent quantities [10], which are constants of motion and which are in involution. 
In fact define the gauge-invariant quantities [11] 



N 



x U ak ' ai 

Then, we can easily show that the 's are conserved quantities, = 0, and that 
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they are in involution 



{F k , Fi} P = 0, (M = 1,..,JV). 



(32) 



N — 1 of them are independent, because 

N 

Fk = x 2 = const.. (5 

fe=i 

The Hamiltonian for the Neumann system is expressed in terms of the F k 's as 



3. The Jacobi system 

The Neumann system is closely related to the mechanical system of the geodesic 
motion of a particle on an ellipsoid, which was already shown to be integrable by 
Jacobi in 1838 [2]. The Lagrangian of the geodesic motion is written as C = 
Y^k=i l'kPk — 7~Li where the Hamiltonian H of the geodesic motion is given by 




(34) 



k=i 




(35) 



with Hj = = \ ^2k=i Pk ■ T ne P r i me " ' " denotes the differentiation with 
respect to the parameter s, e.g. q' k = dqk/ds. The Hamilton equations derived 
from (35) are 




(36) 



(37) 
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N 

QkPk = 



It follows from (36), (38) and (39) that 



N , N N o 

V- Qkq k _ \- qkPk , % _ n 

fe=i fe fc=i fc fc=i fe 



so that 

r? = 0. 

Also we have from (36), (37) and (39) 



N i I N o Wo 

\- 1kP k +P k q k _^P% tST q k - n 

fc=i K fe=i K fe=i fc 



Thus we have 

N „2 



i / i V y \ 



where 

iV 

i? 2 - ]>» 2 ). 



The Dirac brackets are found to be 



{Qi, Qj}D = 0, 

r i r 1 QiQj 

Pj}D = k 3 ~ 

r \ 1 (QiPj ~ gjPO 
te, Pib - ~ W ^ 
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In order to see the connection between the Neumann and the Jacobi systems, 
let us introduce new variables xi, (i — 1, N) and the time variable t by [7] 



** - s|- (48) 
ft - *#■ <«) 

Here k 2 = (R 2 Yld=iPi Ih) 1 can be shown to be a constant. Therefore, we have 
the relation 

N ^ 

s 2 = 1, (50) 



1 ( Pl\ 



R 2 

where s = ds/dt. The nonlinear mapping from q to x by (48) is a kind of the Gauss 
mapping, which maps the ellipsoid onto the sphere with the radius r. 

Differentiating both sides of (48) with respect to t twice, we obtain 



Also differentiating both sides of (49) with respect to t one obtains 



(51) 



S = 2 . (52) 



Substitution of (49), (50) and (52) into (51) gives 

/ N n 2 \- 1 



Ok 



N 2 
x k x k -2 X l \ 

i=i 1 



(53) 
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Putting l/b k = a k an d f 2 — x 2 we finally arrive at Neumann's equations again: 



N 

x k = -a k x k - ^^2(xf - aixf). (54) 

OC 

1=1 



It should be noted that the first derivative of x k can be written as 

/CRN 

^fc = J/fc 2~x k , (55) 

where 

T ( p k T 

Vk = j^j^ s + fix k = — (p k s + fiq k ) (56) 

with \x an arbitrary gauge parameter, which may be a function of x and y. Eq.(55) 
is corresponding to (20). Differentiating (56) with respect to t, one finds 



1 (x ■ y) (x ■ y) 

Vk = —r x k + — —Vk ~ 2/i — — x k + fj, x k + nx k (57) 

Ofe Ob Ob 

This equation agrees with (22), if one fixes the gauge in (22) as follows: 



A = 4-2/!^U/i 2 +A (58) 

Ob Ob 

with a k = l/b k . Needless to say that Neumann's equations also follow from (55) 
and (57). 

It is to be noted that the Jacobi system is also Liouville integrable, because 
there exist iV gauge-invariant quantities G k , (k = l,...,N), corresponding to F k 
in (31) [11], 

with 

Km = q k pi - qipk, (60) 
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which are also shown to be constants of motion, Gj/ = 0, and 



{G k , Gi} P = 0. (61) 



N — 1 of them are independent, because 

iV„ / N 9\ N 2 / N \2 

Ef =(i-£g)£^- E^) -o. w 

fe=l K v k=l kJ k=l k v fc=l k 7 
Hamiltonian for the Jacobi system can be written in terms of G k as follows: 

1 N 1 

^=2^^ = / ( 63 ) 
fc=i 

One finds from (31), (48), (56) and (59) the relation 

,^2 a2 , y . 2 ... 2 ZV^ 

*=^«-*>=-ljr£A' (64) 

from which one obtains useful formulas such as 

N N F, r 2 

fc=l k=l 
N N F, 

£fe = £^ = °- (66) 

k=l k=l k 

It also follows from (50) and (56) that 

N 2 2 

ES-^ + f)' «*> 

It is interesting to note that this equation is transformed to 

N .2 , (,, , ,,\2 



Ef^O^) m 



k=l ak 



under the gauge transformation Zk = Uk + vx k . In consequence, the phase 

space of the Jacobi system is mapped by (48) and (56) onto the sphere defined by 
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x 2 = r 2 and the manifold given by (67) in the phase space of the Neumann system. 
If one fixes the gauge parameter /x, the mapping is obviously bijective. 



4. The Rosochatius system 

In this section we generalize the formalism developed in preceding sections and 
consider the Rosochatius system [9,10] in the generalized framework. Let us start 
with rewriting (48) in the form 

Xi = (71) 

where fi(q), (i — 1, N), are regular functions of the coordinates qj, (j — 1, N) 
which correspond to those introduced in the Jacobi system. The function R = R(q) 
is assumed to be given in terms of fi(q) as follows: 

N 



r 2 = J2^ ( 72 ) 



1=1 



so that 



N o N 

R 2 

i=i i=i 



showing that Xi, (i — 1, N) are variables on the N — 1 sphere S 1 . 
Differentiating X{ with respect to t, one obtains 



Xi = yi y-Xi, (74) 

x z 



where 



T 

Vi = ^ /• + fiXi (75) 



with /i the gauge parameter, s = ds/dt and f[ = dfi/ds. 
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Differentiating both sides of (74) once again with respect to t, one obtains 



As we noticed in the case of the Neumann system, yi is written in terms of X{ and 
yi as follows: 

iji = AiXi + By i: (77) 



where 



If we take 



R = <M 

x 2 ' (78) 
Ai = -(H, (i = l,...,N). 



Ai = -a { + ^ (79) 

x\ 



in addition to B — (xy)/x 2 = 0, then (76) turns out to be 



N N 



x i = -a i x i + ^-^(x 2 -^^1 + ( 8 °) 



x i X k=1 k=l X k' 



which are the equations of motion of a particle on the sphere S N 1 under the 
influence of the potential 

1 N 

( 81 ) 

k=l X ^ 

Eqs.(80) are known as the Rosochatius equations, and the Hamiltonian system 
governed by the potential U(x) given above is called the Rosochatius system [9,10]. 
This system has been shown to be completely integrable by Rosochatius. 
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The geodesic equations (43) in the Jacobi system is rewritten in terms of the 

fi's as 

J= i j 

Our final task is to find the geodesic equations for the Rosochatius system corre- 
sponding to (82) in the Jacobi system. In order to do this we first consider the 
2iV-dimensional Neumann system given by the Hamiltonian 

j 2N 2N 
i=l j=l 

where 

1 N 

7-L2N = g Y^\- W k + W k+N + a k4 + a k+Nzl+N\ ( 84 ) 
k=l 

with Wi = Zi (i — 1, 2iV) and a k+ ^ = a k [3]. 

We introduce new variables x k and 9 k , (k = 1, N), by 

z k = x k cos9 k , z k+ N = x k sin9 k (85) 

and rewrite TC2N in terms of the new variables as 

1 N 

n 2N = - + x ¥l + a k x 2 k ). (86) 

k=l 

We then restrict ourselves to the specific case in which the angular momenta 
z k w k+ N — Zk+N w k> (k = 1, N) are the constants of motion, that is, 

ZkWk+N - z k+N w k = x\b k = ^ (k = 1, N) (87) 

with Cfe, (k = 1, ...,JV) real constants. Substituting (87) into (86), we obtain the 
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Hamiltonian for the Rosochatius system 

1 N 

^^E^ + V fl ^4 (88) 

Z k=l X k 

We next consider the 2iV-dimensional Jacobi system. The Hamiltonian is given 
by (35) with N replaced by 2N and b^ = &fc+jv = 1/ajfc, (k = 1, iV). The geodesic 
equations are then written as 

where 



1 -r? 



2 i = ^j. ^-+jv = ^^+iv- ( 90 ) 
In view of (85) it is convenient to introduce the variables fa, (i — 1, N) by 

ft = fa cos6»i, ft +iV = fa sm9i. (91) 

Substituting (91) into (89), we obtain 

1 N 1 

- m 2 + J2 -t(/j) 2 + /M) 2 ] = °> ( 92 ) 



VlO'i + = 0. (93) 
Integrating (93) we find that ffO^ are constants of motion. We set 



fM = ^r- 04) 



and if fj = qj/bj as in the iV-dimensional Jacobi system, we finally obtain 
d 2 qi di 1 ^> 1 \fdqj\ 2 i dj 



ds 2 qf R 2 <-[0jL\ass qp 



with R 2 given by (44). 
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It is interesting to note that N integrals of motion H k , (k = 1, ...,N), for the 
Rosochatius system are given by 

Hk = Fk + Fk+N = 4 + 4 +N + L 2 ^Ji^ + l- 2 J2 ^rf^' (96) 

where Jfcz = z k wi—ziw k . It should be noted that both F k and F^ + tv contain singular 
terms because the Hamiltonian (83) is degenerate, a k = a k+ N,(k = 1,...,N). 
However, such singular terms cancel in F k + F k+ n and we have 

N 



with[9,10] 



x U ak ~ ai 



2 2 

Hi = {x m - x iyk f + ^ + (98) 

In a similar way we obtain N integrals I k) (A; = 1, N) for the equations (95) 
combining pairs of integrals for the 2iV-dimensional degenerate Jacobi system 

2iV V 2 2iV K 2 

5£ - 6/ ^ b k+N - h 

where ii" fcZ = ( k 7Ci - (iir k with 

Ck = h9k = b k f k cos9 k = q k cos6 k , 

Ck+N = b k+ N9k+N = ?fc sm ^fc) 

7T k = C' k = (p k cos6 k - q k 9' k sme k ), 
Kk+N = Ck+N = (Pksin6 k + q k 9' k cos9 k ). 

Substituting (100) into (99), we find 

(101) 



(100) 
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with 

.2 j Jl 



4 = («M-«P*) 2 + % + % (102) 

Ik 1l 



Since it can be shown that 



^ = 0, (k = l,...,N), (103) 



Ik are indeed conserved quantities of the system governed by (95). 



5. Summary and concluding remarks 

In the first part of this work we dealt with the Neumann and the Jacobi systems 
in the classical framework of constrained Hamiltonian systems. We calculated 
Dirac brackets for canonical variables. We noticed that a transmutation from the 
second-class constraints to the first-class ones occurred by changing dynamical 
variables appropriately. As a consequence, both systems acquired gauge freedom 
in terms of residual gauge-parameters. 

We focused on the relationship between the Neumann and the Jacobi sys- 
tems. The mapping from the phase space of the Jacobi system to that of the 
Neumann system was executed by the Gauss map. Affine connections appearing 
in the geodesic equations in the Jacobi system are given by 

T * = wlfc- < 104 > 

However, these affine connections do not satisfy the identities which genuine affine 
connections should satisfy, because the geodesic equations are not all independent. 

In the second part of the work we considered the Rosochatius system. We 
derived the Hamiltonian 7{r for the Rosochatius system from the 2A^-dimensional 
Neumann Hamiltonian. The geodesic equations for the Rosochatius system should 
be obtained by making use of the 2A^-dimensional Jacobi equations. However, the 
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final result given by (95) is by no means the geodesic equation for qi, because (95) 
does not have the standard form of the geodesic equation, which is usually written 

as 

d 2 fk r k dfi dfj = Q 
ds 2 %J ds ds 

N (105) 

2^1=1 fa ij=l 

where is the regular function of fi, (i = 1, ...,JV), fa = d(p/dfi and faj = 
d^/dfidfj. 

However, we have 



d 2 6 k , 2 % d6> fc 



+ 



qi~ ds ds 
(d9 k \ 2 



0, 



N 



1 ^ 1 

3=1 J 



l\ ds 



\ds) . 



(106) 



from (93), (94) and (95). These equations are geodesic equations in the 2N- 
dimensional space described by the coordinates (<%, 9k), (k = 1, ...,N). 

We found that (95) are integrable and (k — 1, N) are conserved quanti- 
ties of the system governed by (95). In other words, the set of equations given by 
(95) constitutes an integrable system dual to the Rosochatius system. 

There are some important problems on the Rosochatius system with its dual 
system left unsolved, such as finding classical solutions in arbitrary dimensions and 
the problem of quantizing the Jacobi-Neumann and the Rosochatius systems. Some 
classical solutions to the Neumann system have been found by several authors [6]. 
Quantization of the Neumann system is discussed in some extent by Gurarie [14]. 
We would like to discuss extensively those remaining problems in a forthcoming 
paper. 
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